Coupling between scattering channels with SUSY 
transformations for equal thresholds 

December 23, 2008 

Andrey M Pupasovj: 12 , Boris F Samsonov 1 , Jean-Marc 
Sparenberg 2 and Daniel Baye 2 

1 Physics Department, Tomsk State University, 36 Lenin Avenue, 634050 
Tomsk, Russia 

2 Physique Quantique, CP. 229, Universite Libre de Bruxelles, B 1050 
Bruxclles, Belgium 

E-mail: pupasov@phys.tsu.ru, samsonov@phys.tsu.ru, jmspar@ulb.ac.be, 
dbaye@ulb . ac . be 

Abstract. Supersymmetric (SUSY) transformations of the multi-channel 
Schrodinger equation with equal thresholds and arbitrary partial waves in all 
channels are studied. The structures of the transformation function and the 
superpotential are analyzed. Relations between Jost and scattering matrices of 
superpartner potentials are obtained. In particular, we show that a special type 
of SUSY transformation allows us to introduce a coupling between scattering 
channels starting from a potential with an uncoupled scattering matrix. The 
possibility for this coupling to be trivial is discussed. We show that the 
transformation introduces bound and virtual states with a definite degeneracy at 
the factorization energy. A detailed study of the potential and scattering matrices 
is given for the 2x2 case. The possibility of inverting coupled-channel scattering 
data by such a SUSY transformation is demonstrated by several examples (s — s, 
s — p and s — d partial waves) . 



1. Introduction 

The inverse problem in quantum scattering theory plays an essential role in theoretical 
and mathematical physics. Following the needs of soliton theory [1, 2] and elastic 
scattering processes in atomic and nuclear collisions [3, 4], the one-channel inverse 
problem reached a high perfection level of development summarized in a number 
of books (see e.g. [5, 6]). The most promising way of solving this problem from a 
pragmatic viewpoint is related with the so called method of Darboux transformations 
[7] which is also known under the name of supersymmetric (SUSY) quantum mechanics 
[8, 9] (for a review of inverse scattering with SUSY quantum mechanics, see [10]). 
In contrast, the multichannel inverse scattering theory, which describes inelastic 
scattering processes [3, 4], only gained a fragmental level of development (see e.g. 
[6, 10]). In particular, after a few papers devoted to an application of SUSY 
transformations to multichannel inversion published by Cannata et al between 1988 
and 1992 [11, 12, 13, 14, 15] and two papers published in 1997 [17] and 2000 [18], 
where phase equivalence is discussed, a publication gap occurred till 2006. We think 
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that the main reason for this inactivity is the lack of a detailed analysis of properties of 
general SUSY transformations applied to the scattering inversion in the multichannel 
case. 

For potentials with different thresholds such an analysis permitted us to introduce 
non-conservative SUSY transformations [19], with the help of which, starting from a 
decoupled potential, new exactly solvable coupled potentials were obtained [20, 21, 22]. 
We were also able to give a careful analysis of the scattering properties of the known 
Cox potential [23] and apply it to construct a realistic model of the magnetic Feshbach 
resonance phenomenon [24]. For equal thresholds, to the best of our knowledge, only 
two papers exist [11, 14], where only a particular type of SUSY transformation is 
studied, i.e. a transformation removing a bound state. 

In the current paper, with the equal threshold inverse problem in mind, we rather 
concentrate on the necessary ingredients for a single SUSY transformation to generate 
coupled scattering matrices, starting from a decoupled potential. We first show that 
such a transformation introduces bound and virtual states at the same energy and we 
calculate their degeneracy. Next, we discuss the possibility to get a trivially or non 
trivially coupled scattering matrix when both the potential and Jost matrix are non 
trivially coupled. By trivially coupled scattering and Jost matrices, we mean non- 
diagonal matrices which may be diagonalized by energy-independent transformations. 
Similarly, a multichannel potential is nontrivially coupled if its matrix cannot be 
diagonalized by an r-indcpcndent transformation, where r is the radial coordinate. 
Below, we will be able to answer the following questions: does a non trivial coupling 
of the potential imply a non trivial coupling of the scattering matrix? Does a non 
trivial coupling of the Jost matrix imply a non trivial coupling of the scattering matrix? 

The paper is organized as follows. In the next section, we fix our notations 
regarding SUSY transformations as applied to multichannel scattering theory. In 
Section 3, we analyze the asymptotic behaviour of the superpotential for SUSY 
transformed potentials and establish conditions to get potentials with non-trivially 
coupled 5-matrices from uncoupled potentials. The analysis of the Jost-matrix 
determinant permits us to find modifications produced in the spectrum by SUSY 
transformations. The rest of the paper is devoted to a detailed analysis of the most 
important case from a practical point of view, i.e. the case of two-channel systems 
which can model a rich class of physical systems [25]. In Section 4, the long range 
behaviour of the transformed potential is studied. In particular, we indicate an unusual 
possibility for a SUSY transformation consisting in the interchange of the centrifugal 
terms by the scattering channels. In section 5, we find the Jost and scattering matrices, 
phase shifts and mixing angle for the SUSY partner potential. Section 6 first presents 
examples of coupled potentials with trivially coupled scattering matrices. Then we 
exemplify coupled potentials (s — s, s — p and s — d partial waves) with non-trivially 
coupled scattering matrices and analyze their eigenphase shifts. In the Conclusion, a 
summary of the obtained results is given. 

2. Outline of scattering theory and SUSY quantum mechanics 

Our starting point is the N x N multi-channel Schrodingcr equation (E = k 2 ) 
HiP(k,r) = k 2 ^(k,r), H = -I N f^ + V(r), V = + V(r) , (1) 
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where r is the radial coordinate, r G (0, oo), V(r) is an N x N real short ranged 
symmetric potential matrix, I — diag(7 l7 . . . , l^), lj = 0, 1, . . ., is a diagonal matrix 
of angular momentum quantum numbers, and ip may be either a matrix- valued or a 
vector- valued function. Here and in what follows we also denote Z±l (and v±l below) 
matrix I ± 7jv {y ± -fjv) where /at is the N x N identity matrix. Matrix t/ determines 
the singularity of the potential at the origin 

V(r ->0) =r- 2 v{v + l) + 0{l) 1 v = diag [v u v 2 , . . . , i>jv] , ^ = 0, 1, • • • • (2) 

Note that V does not contain a Coulomb-like r _1 singularity. Matrix I defines the 
asymptotic behaviour of the potential at large distances 

V(r -> oo) = r- 2 l{l + l)+o{r- 2 ) (3) 

which is typical of coupled channels involving various partial waves. 

We define the regular N x N matrix solution ip(k,r) of (1) according to [3] 

<p(k, r -» 0) - diag(r" 1+ 7(2i/i + 1)!!, . . . , r v " +1 /(2v N + 1)!!) 

=: r' /+1 [(2t/ + l)!!]- 1 . (4) 

The Jost solution f(k,r) has the usual exponential asymptotic behaviour at large 
distances 

f(k, r oo) -> diag[e 4fcir , . . . , e lfclvr ] := e lfcr . (5) 

The Jost solutions f(k, r) and f(—k, r) form a basis in the matrix solution space. The 
regular solution is expressed in terms of the Jost solutions as 

<P(k, r) = ^ [f(-k, r)F(k) - f(k, r)F(-k)} , (6) 

where F(k) is known as the Jost matrix. From (4) and (6) it follows that 

F(k) = lim [/ T (fc, r)r v ] [{2v - l)!!]" 1 . (7) 

r — >0 

The physical solution, which appears in the partial wave decomposition of the 
stationary scattering state, behaves as 

*(fc, r -> oo) oc e il %~ ikr - e - a ^ +lkr S{k) . (8) 

The scattering matrix S(k) is expressed in terms of the Jost matrix as 

S(k) = e a %F{-k)F- 1 {k)e a % . (9) 

Below we will consider SUSY partners V(r) of the potential V(r) and find the 
corresponding partners for both the Jost and scattering matrices. 

For the reader's convenience, we now describe the general scheme of SUSY 
transformations. It is known that solutions of the initial Schrodinger equation (1) may 
be mapped into solutions of the transformed equation with help of the differential- 
matrix operator L, 



tp(k, r) — Lip(k, r) :- 



w(r) - In4~ 
dr 



1>(k,r). (10) 



The transformed Schrodinger equation has form (1) with a new potential 

V(r) = V(r)-2w'(r), (11) 

Here the prime denotes the derivative with respect to r. Matrix w(r), called the 
superpotential, is expressed in terms of a matrix solution u of the initial Schrodinger 
equation 

Hu(r) = -K 2 u(r) , (12) 
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as follows 

w(r) =u'(r)u- 1 (r). (13) 

Here n is known as the factorization wave number. It corresponds to an energy 
E = —k 2 , called the factorization energy. Solution u(r) is called the factorization 
solution or (matrix-valued) transformation function. 

Let us also mention that, since Lu{r) = 0, solution r/>(i/c, r) of the transformed 
equation corresponding to the energy E = — k 2 is found as matrix <E>(r) = 
Moreover, this matrix, when chosen as the transformation function for the next 
transformation step, cancels the previously produced potential difference. This means 
that it corresponds to a transformation in the opposite direction. 

In the next section we introduce a transformation which mixes scattering channels 
in an uncoupled potential 

V(r) = V d (r) := diag (V d(1 ,V dfl , V d , N ) (14) 

with given values of v and I. Such potentials may be obtained, for instance, by proper 
chains of usual 1-channel transformations of the zero initial potential (see e.g. [26]) 
resulting in generalized Bargmann potentials (for more details see e.g. [27]). 



3. Coupling transformations 

In the most general case the transformation function may be expressed in terms of 
the Jost solutions as follows 

u c (r) = fd(-iK, r)C + f d (i K , r)D , (15) 

where the real constant matrices C and D are arbitrary. Here and in what follows 
subscripts d and c stand for quantities related with diagonal (uncoupled) and non- 
diagonal (coupled) matrices, respectively. To obtain a Hcrmitian potential after a 
transformation with transformation function (15) matrix C T D should be symmetric 
[20] 

D T C — C T D . (16) 
We need the asymptotic behaviour of the superpotential w coo := lim w c (r) to find 

the transformed Jost solution and, hence, the Jost and scattering matrices. As it was 
shown in [20] for different thresholds, this behaviour of the superpotential depends 
crucially on matrix C. Below we shortly discuss the method developed in [20] while 
making necessary changes to adjust it for the case of equal thresholds. 

Matrices C and D with a maximal number of independent parameters 
guaranteeing the Hermitian character of the superpotential (13) have the following 
canonical form, 

where Xq = Xq is a real symmetric nonsingular M x M matrix, and Q is an 
(N - M) x M real matrix so that rankC = M. 

The asymptotic matrix w coc is determined by the behaviour of transformation 
function (15) at large distances 

Mr -™^ A { I N _ u e-*r ) ' \ Q In-m)- (18) 
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From (13) and (18) we obtain 

—In-m 



w coo = kA[ ° U" 1 (19) 



with 



Tf. 

(20) 



4 -i = a t Im + Q'Q 

V In-m + QQ t 

Comparing this result with that obtained in [20], we conclude that the main 
difference between equal and different thresholds is the non-diagonal character of the 
superpotential at infinity. Note that superpotential w coo has a richer structure than 
that previously reported by Amado et al [11]. Their result corresponds to the choice 
M = 1 when w coo is expressed in terms of a single (N — l)-vector Q = (qi, . . . , qN-i) T ■ 

Once w coo is determined one can calculate the Jost solution f c (k, r) and the Jost 
matrix F c (k) for the transformed potential V c . The Jost solution takes the form 

f c {k,r)=L c f d (k,r){w coo -ikI N y 1 , L c := w c (r) - I N -^- . (21) 

ar 

The factor (w coo — ikl^)^ 1 is introduced to guarantee the correct asymptotic behaviour 
of f c (k,r) (see (5)). In order to find the Jost matrix, we first consider the behaviour 
of the superpotential in a vicinity of r — which depends on the character of the 
transformation solution (15). Below we will assume that there is no bound state 
at the factorization energy, detF d (in) ^ 0, and each column of the transformation 
solution is singular at the origin. Using the property 

/d(-fc,r - 0) = f d (k,r^ 0)F d (-k)F- 1 (k) + o(rn (22) 

which follows from (6) and the invertibility of F d (ihi), one finds the behaviour of the 
transformation solution at the origin, 

u c (r -> 0) = f d (iK, r)[F d {-iK)F^{iK)C + D] + o{r v ) . (23) 

We assume 

Aet(F d {-in)F^ 1 {m)C + D) ^0, (24) 

which can always be provided by a proper choice of matrices C and D. The leading 
term of the superpotential at r — > reads 

w c (r -> 0) = -r~ x v + o(l) , vj>0, (25) 

where we used the Laurent series for the Jost solution 

f(k, r) = r- v {2v - l)!!F T (fc) + r-" +1 &i(fc) + o(r-" +1 ) . (26) 

It follows from (2) and the Schrodingcr equation that b\(k) — 0. We note the diagonal 
character of superpotential (25) at the origin. The singularity at the origin of the 
transformed potential, 

V c (r - 0) -> r- 2 D{v + 1) = (V d - 2w' c )\ r ^ = r" 2 ^ - 1), (27) 

decreases by one unit, v — ► v = v— 1. Hence we can apply our coupling transformation 
to potentials for which matrix v is positive definite, v > 0, a property we will assume 
to hold throughout the paper. 

The Jost matrix can be obtained from expression (6) of the regular solution 
(p c {k, r) corresponding to V c . The regular solution of the transformed potential <p c (k, r) 
is determined by (4) with the singularity parameter v. To derive it, we act on 
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both sides of expression (6) of the regular solution ip d (k,r) for potential V d with 
the transformation operator L c . From (4), (10) and (25), it follows that 

L c tp d (k, r) = -<p c (k, r). (28) 

Taking into account (21), we rewrite (28) as 

<p c (k,r) = -— [f c {-k,r){w coo +ikI N )F d (k) - f c {k,r)(w coo - ikI N )F d {-k)] . (29) 

Comparing (6) and (29) we find a relation between the initial and transformed Jost 
matrices 

F c (k) = -(ikI N + w coo )F d {k) . (30) 

For M — and M = N, Q is absent and A = I N in (19). When M = N, 
the superpotential at infinity (19) becomes proportional to the identity matrix, 
w coo — kIn. The transformed Jost matrix (30) becomes diagonal. Similarly, the 
case M — leads to w coo = —kIn- From here we draw an important conclusion. 
The necessary (but not sufficient) condition for a non-trivial coupling in the Jost and 
hence scattering matrices is < M < N. This will be assumed in the following. 

As already mentioned, a non-trivial coupling in the Jost matrix requires not only 
a non-diagonal Jost matrix, but also the impossibility to diagonalize this matrix by 
a fc-independent transformation. It is clear that matrix ikljq + w coo from (30) can 
be diagonalized by a fc-independent transformation. When F d is not proportional to 
the identity matrix, channels in F c (k) are coupled in a non-trivial way. Nevertheless, 
this property does not guarantee the non-triviality of the S'-matrix. As it follows from 
definition (9), the 5-matrix will be trivially coupled when the product F d (—k)Fl L 1 {k) 
is proportional to the identity matrix, i.e. when 

F d ;j(k) = \F d .j(k)\e- i5 ^ S d , 3 {k) = {-l)he»sW , 

j = l,...,N. (31) 

In particular, if all single-channel potentials have the same S-matrix, i.e. if they 
are phase equivalent [28] or isophasc [29], the S- matrix resulting from a SUSY 
transformation keeps being trivially coupled. When the initial 5-matrix is not 
proportional to the identity matrix, one may expect non-trivially coupled channels. 

The analytic expression (30) for the Jost matrix allows us to study spectral 
properties of the transformed potential (11) [21]. The positions of the bound/virtual 
states and resonances arc defined as the solutions of detF c (fc) = 0. As it follows from 
(30) and (19), the Jost-matrix determinant is given by 

detF c (fc) = (-l) N (ik + K ) M (ik-K) N - M detF d (k) (32) 

since w coo has the M fold degenerate eigenvalue n and the N — M fold degenerate 
eigenvalue — k. Therefore, if det F d (k) has no pole at k = ±in (this property is 
assumed to hold in the rest of the paper), the SUSY transformation leads to a new M 
fold degenerate bound state with kb — in, = —n 2 and an N — M fold degenerate 
virtual state with k v = —in, E v = —n 2 . 

Now we continue to compare our method with the approach developed by Amado 
et al [11]. For that, we calculate the asymptotic behaviour of matrix = [u c (r)^]~ 1 , 
which upon using (18) reads 
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The M first columns of $(r) are vectors decreasing at infinity. According to 
(23) and (26), $(r) is a regular solution, $(0) = 0. Therefore these vectors 
correspond to the bound state wave functions of the coupled system appeared 
after the SUSY transformation. This confirms that the energy level of this bound 
state is M fold degenerate. All the other columns in $(r) correspond to virtual 
states. For the particular case M — 1 this asymptotic form just corresponds to the 
transformation function used by the authors of [11] for decoupling a coupled problem. 
We thus conclude that their transformation corresponds to a particular case of our 
transformation when realized in the opposite direction. 

Another useful remark is that although the superpotential w c (r) depends on 
parameters Xq, the Jost matrix F c {k) and, hence, the S-matrix are Yo-independent. 
This means that superpotential w c (r) leads to a family of potentials, parameterized 
by the entries of Y , having the same scattering properties. 

Below we concentrate on the two-channel case with equal thresholds and arbitrary 
partial waves. The coupling SUSY transformation produces in this case one bound 
state and one virtual state. First we will analyze the long range behaviour of the 
transformed potential. 



4. Long range behaviour of the transformed potential 

In the two-channel case, according to (3), the initial diagonal potential has the 
following long-range behaviour 

v *<r^<*>)^-?\ o Hh + i))- (34) 

The Jost solution at large distances is expressed in terms of third kind Bessel functions 
Hj^^z), also called first Hankel functions (see [30] for a definition) 

/ d (fc,r-c») = diag[MAr),Mfcr)] . hi{z) = i l+1 {nz/2)i (z) . (35) 
The recurrence relations for hi(z) and its asymptotic behaviour 



(36) 



follow from those for h[ 1] ' (z) [30]. 

For the coupling transformation, according to our previous discussion, we choose 
transformation function (15) with matrices 

which contain only two independent parameters x and q. The restriction on the 
parameters 

x + F dA (-iK)F£(iK) - q 2 F di2 (-iK)F^{iK) + (38) 
follows from (24). Transformation solution (15) reads 

„ ( r \ _ f fd:i{-iKr)+xf d a{iKr) -qf d -i{inr) \ , . 

Let us consider the first two terms in the asymptotic behaviour of the 
superpotential (13), w c (r — > oo) = w coo + W-ir -1 + o(r _1 ). The first term w coo 



Coupling between scattering channels with SUSY transformations 



8 



has been calculated for an arbitrary number of channels in Section 3. Thus from (19) 
we obtain 

— if? ( 1 /! 1 ) ■ <«> 

Another parametrization for w coo is useful, 

/ cos a sin a \ a 
w cod = k\ . , g = tan — . (41) 

co ° ^ sin a - cos a J ' y 2 

Note that a non-zero value of W-\ will lead to a modification of the long range 
behaviour of potential (11) with w' c (r — > oo) = — w_ir~ 2 + o(r~ 2 ). 

In order to establish the asymptotic behaviour of the potential V c (r) = Vd(r) — 
2w' c (r), wc replace fdj(iiKr) in (39) by its asymptotic form given in (35) and neglect 
in (11) and (13) all exponentially decreasing terms such as /i^ (mr)/ii 2 (mr). Taking 
into account 

/ ^^1 (^Kr)/^ ^2 (-^Kr) = l + ^-[/ 1 (^ 1 + l)-/ 2 (/ 2 + l)]+o(r- 1 ), r oo, (42) 
combining (11) with (34) and using parametrization (41), one finally gets 

1 ( h(h + ±) 



v °<r^°°)-^\ o i 2 {h + i) 

[^2(^2 + 1) - h(h + l)]sina / sina -cosa\ 

r 2 \ — cos a — sin a J ' 

A similar asymptotic behaviour of the matrix potential is obtained from the Gelfand- 
Levitan equation in [31]. 

From (43) we conclude that, for Zi 7^ I2, the transformed potential has a non- 
zero long range coupling, V c; i2 ^ 0. Moreover, it is impossible to associate diagonal 
entries of V c given in (43) with usual centrifugal terms. The only way to avoid this 
inconvenience is to fix a = ±ir/2 (or equivalently q = ±1), which leads to a physically 
reasonable long range behaviour 

t -i,-£(«v 11 («) 

Having compared (34) and (44) we find the modification of the corresponding angular 
momentum quantum numbers under the SUSY transformation I = diag(Zi, I2) — > I = 
diag(?2, h)- For short, this unusual property of the SUSY transformation will be called 
the exchange of the channels angular momenta. Summarizing, we get an additional 
constraint q — 1 (the dual case q = — 1 leads to the same transformed potential except 
V c -i2 — * —V c -\2) to consider only physically reasonable potentials in the case I2 ^ h- 

5. Transformed Jost and scattering matrices, eigenphase shifts and 
mixing angle 

In the two-channels case, introducing w coc as given in (41) into (30) provides an 
explicit relation between the transformed Jost matrix F c (k) and the initial diagonal 
Jost matrix Fd(k), 

F c (k) = -( lk + KC ° Sa ., KSina W). (45) 
w V « sin a ik — k cos a J 
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From (32), we obtain det F c = (fc 2 + K 2 )det F d . The coupling transformation produces 
one bound state and one virtual state, in agreement with the general properties of the 
transformed Jost matrix analyzed in Section 3. 

Once the transformed Jost matrix F c (k) (45) is found, the S- matrix may be 
obtained according to its definition (9), where we have to take into account the change 
of attribution of the angular momenta Z — > Z by the coupling transformation, 

S c (k) = {-ikh +w coo ){-l) l S d {k)(ikI 2 +w coo )- 1 e il % . (46) 

The diagonal matrix 

S d (k) = di&g{e 2l5d ^ k \e 2lSd ^) (47) 

is obtained from the diagonal Jost matrix F d (k) before the transformation. One can see 
that for the particular case of identical partial waves, I cx I 2 , our result (46) reproduces 
the corresponding relation (17a) from [11]. For different partial waves however, the 
modification of the angular momenta leads to the appearance of additional phase 
factors e*'5 and (—1)'. 

Let us now find the transformed cigenphase shifts 5 c .j(k), j ~ 1,2, and the mixing 
angle e(fc). Since the scattering matrix is unitary and symmetric, there exists an 
orthogonal transformation R(k) diagonalizing this matrix, 

Rc {k) { S ca2 (k) S c , 22 {k) ) Rc ^ ~ (, e 2 ^« ) -( 48 ) 
The rotation matrix R c (k) is parameterized by a mixing angle e = e(fc), 

^-(-"18) ctS!)- ,49) 

which is expressed in terms of S'-matrix entries as 

tan2 £ (fc)= gpH^j (50) 

Here one can distinguish three essentially different cases: 

(i) the difference between the angular momenta is odd, l 2 = (li + 1) (mod 2) ; 

(ii) the difference between the angular momenta is even, l 2 ^ h, l 2 = l\ (mod 2); 

(iii) the angular momenta coincide, l 2 = h- 

Note that case (i) does not correspond to any reduction of the rotationally 
invariant three-dimensional scattering problem, since in this case any nontrivial 
coupling means a parity breakdown (see e.g. [4]). For the sake of completeness we 
will analyze this case also although the corresponding system of coupled equations 
has no direct relation to a scattering problem. Moreover, we will use the usual 
scattering theory terminology in this case also, although from the point of view of 
a three-dimensional scattering it bears only a formal character. In cases (i) and (ii), 
we should put a — tt/2, which is not necessary in the third case (see (43)). 

Definition (47) allows writing 

S d -i(k) + S d . 2 (k) = e 2«*i(*) +e 2«d i2 (fc) = 2e^ 1+ ^ 2 >cosA, A = S d . 2 -S d:1 ,(51) 

and (50) leads to expressions for the mixing angle in the three cases as 

„ ,,. 2(-l)(' 2 ~ Zl ~ 1 )/ 2 Ksina:(fcsin A — KCOsacosA) 

i : tan2e fc = -i '-— v — '- , 52 

2nk cos a cos A - (fc 2 - k 2 ) sin A 
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.... ..... .,. 2(-l)(' 2 ~ (l )/ 2 Ksina («cosasin A + fccos A) . . 

(11), (111) : tan2e(fc) = : — — — — — - . (53) 

v / v ' w sinA(fc 2 - K 2 cos2a) - 2KfccosacosA v ; 

Since q = 1 and a = n/2 in cases (i) and (ii), expressions (52) and (53) are simplified 
to 

(i) : e(fc) = (-l)(^-«i-i)/2 arctan * ; ( 54) 

K 

(ii) : tan2e(fc) = (-l)C 2 -' 1 )/ 2 cot A. (55) 

k + k 

We will assume below that the scattering matrix of the transformed potential 
satisfies the effective range expansion (see e.g. [32]), which implies 

cot 5 c;h2 (k -> 0) = a h2 k'^ +1 ^ + o (fc-^^+D) , 

e(k -» 0) = eofc 1 ' 2- ' 1 ' + o (fc 1 ' 2- ' 11 ) . (56) 

Since there are rather simple analytical expressions for the mixing angle, we will 
analyze restrictions on parameters of the SUSY transformation which follow from the 
second equation in (56). 

In case (i), (54) satisfies the effective range expansion (56) when \l 2 — h\ = 1 and 
violates (56) when|Z2~^i| > 1. The important property of the coupling transformation 
in case (i) is that the transformed phase shifts coincide with the initial phase shifts, 
i.e., 

R T c {k)S c (k)R c {k) = Q e2 Z i(k) J . (57) 

Therefore, one may separately fit the phase shifts for the l\ and l 2 waves before the 
coupling transformation. 

In case (ii), the effective range expansion for mixing angle (55) leads to the 
restriction cot A = or <5d.2 (0) — f>d,i(0) = (n + l/2)7r . According to the Lcvinson 
theorem (see e.g. [3]) this means that the potential Vd supports a bound state at zero 
energy. 

Finally, in case (iii) there is no any additional constraint since e(fc — > 0) = const. 

Having established properties of the transformed phase shifts and the mixing 
angle, we will consider in the next section some schematic examples of scattering for 
the s — s, s — p and s — d coupled channels. 



6. Examples 



To illustrate the difference between couplings in potential, Jost and scattering 
matrices, we construct in this section nontrivially coupled potentials having trivially 
coupled S'-matriccs and both trivially and non-trivially coupled Jost matrices. After 
that we exemplify SUSY transformations leading to non-trivially coupled S'-matrices. 
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6.1. Coupled potentials with uncoupled S-matrices 



V(r;j3) = — 2-j-2 In W [sinh(K r), sinh(/-c 2 r), exp(Kir) + /3exp(— Kir)] , (58) 



Let us consider the 1-channcl potential expressed in terms of a Wronskian as 
d_ 
dr 2 

k < Ki < k 2 , [3 < -1 , 

which can easily be obtained from the zero potential with the help of the usual (i.e. 
1-channel) SUSY transformations. This potential has one bound state at energy 
E = — k\ and its Jost function has the form 

F(k) = i(k - iKi) \(k + in )(k + iK 2 )] _1 . (59) 

All potentials from the /3-family (58) have the same Jost and scattering matrices. 
Therefore, we can construct a diagonal potential Vd{r) = diag [V(r; fli), V(r; /3 2 )] with 
a two fold degenerate bound state at energy E = —k\. Both its Jost and scattering 
matrices are proportional to the identity matrix 

F d (k) = F{k)I 2 , S d (k) = F{-k)F-\k)h ■ (60) 

As a result, the Jost matrix (30) obtained after the coupling transformation can be 
diagonalized by the same k- independent transformation as the superpotential w coc . 
This just corresponds to a trivial coupling in both Jost and scattering matrices. 

For the coupling transformation we choose the transformation function (39) where 
Jost solutions f d (iK,r) of the Schrodinger equation (26) with potential V d are used. 
To avoid a singularity at finite distance in the transformed potential we impose the 
restriction n > k 2 > «i- Such a transformed potential is shown in figure 1(a). The 
function o~(r) = V r c: i 2 /(V r c;22 — V c; n) demonstrates the non-triviality of the transformed 
potential matrix. If a is a constant, the potential matrix is globally diagonalizable. 
As we see from figure 1 (b), a non constant a means that the transformed potential 
has a non-trivial coupling. At the same time, the mixing angle (53) in the scattering 
matrix is just a constant for A = 



a 



e(k) = --. (61) 



The phase shifts for this potential read 

S c -i(k) = 2n — arctan — + arctan — , (62) 



1 k k 



v-^. k k 
5 C :2(k) = it — 2_, arctan arctan — . (63) 

From here one can see that after the coupling transformation the additional bound 
state increases the value of the phase shifts at zero energy in agreement with the 
Levinson theorem. 

To show the restrictive character of the requirement for the S matrix to be non- 
trivially coupled, we construct below a potential with non-trivially coupled potential 
and Jost matrices but a trivially coupled S'-matrix. This possibility is based on the 
fact that in the single-channel case two different Jost functions may correspond to 
the same scattering matrix [16, 29]. In this case, the two Jost functions differ from 
each other by a real factor for real fc's (see (31)). Therefore if we apply our coupling 
transformation to the following uncoupled system 

V d (r) - diag^r), V 2 (r)] , F d (k) - diag [F^k), F 2 (k)} , S d (k) = S(k)I 2 , (64) 
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Figure 1. (a) Exactly solvable potential V c = V c obtained from two 
uncoupled potentials (58) (/3i = — 2,f3 2 = — 1.5, Kq = 1 , «i = 2.5 , Kz = 
3.5) by the coupling transformation with parameters q = 0.5, k = 6, x = 25. 
(b) Ratio a(r) = V c; i 2 (r)/ (K ;22 (r) - K ; ii(r)). 



from (64), (11), (30) and (46) we can see that the transformed potential and Jost 
matrices cannot be diagonalized by a constant rotation whereas the scattering matrix 
becomes diagonal after the same fc-independent rotation which diagonalizes w coo . 

An example in which we get a non-diagonal Jost matrix and a trivially coupled 
S'-matrix after applying the coupling transformation follows from (64) where we choose 

i(k — ikq) —i 



F d (k) = diag 



(k + ini)(k + i/t 2 ) ' {k + iKo){k + in\){k + 1K2) 



(65) 



S d (k) = ( fc + iK o)( fc + iK i)( fc + iK 2) ^ ^ 
(k — iKo)(k — iKx)(k — in-i) 

Here Ko, k± and are arbitrary real parameters. Matrix v for the corresponding 
potential 

V d (r) = diag [V(r, (3 < -1), V(r, - -1)] , (67) 

is v = diag(l, 3) meaning that v — 1 > and we can apply the coupling transformation. 
Here the non-trivially coupled transformed Jost matrix (30) with as given in (65) 
leads to the following trivially coupled ^-matrix 

n /n ruT ^2 (fc + ««())(£ + i«l)(fc + i/e 2 ) 

S c (k) = (rkh - Wcx ) {k _ lKQ){k _ lKi){k _ lK2){k , + k2) ■ (68) 

The corresponding phase shifts are given by (62) where S c -2 — * 6 c -2 — ft and the mixing 
angle is given by (61). 



6.2. Coupled s — s partial waves 

Using our general scheme described in section 5, we can study the behaviour of the 
phase shifts for the s — s coupled potential. Since the angular momenta in both 
channels coincide, we have here the case (iii) discussed above. Parameter q is not 
fixed from the long-range behaviour of the potential and the mixing angle is given 
by (53). The analysis of this expression is based on the low-energy behaviour of the 
phase shifts dd-,i,2 before the coupling transformation 

e 2iSd ^ 2 = l-2zai, 2 fc + o(fc), (69) 
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where a\ and a 2 are the scattering lengths for each channel. Combining (53) and (69) 
we get 

2(1 + (ai — an) k cos a) sin a , „ , 

tan 2e fc -> = - V 1 2 y r — ^— + o(fc) . 70 

2 cos a + (ai — a 2 )ft cos 2a 

The expansion of the eigenvalues of the transformed scattering matrix at low energies 
reads 

e 2W Ci i,2 = i _ ik ( ai + a 2 ) ± ^(02 - ai) 2 + 4 (l//c - (a 2 - ai) cosa) /k + o(fc).(71) 

An important result from the point of view of inverse scattering corresponds to 
a coupling vanishing at low energies, i.e. when e(0) = 0. This leads to an additional 
link between the parameters, 

cosa =7 — r— , (72) 

(a 2 - ai)K 

where we have used (70) and (56). Hence (71) simplifies into 

e 2**c,i,2 = 1 _ 2ia 2 ,ik + o(k) . (73) 

In this case, the scattering lengths for the transformed potential coincide with the 
initial scattering lengths 01 and a 2 . This property allows us to fit low energy scattering 
data for uncoupled channels thus simplifying essentially the inverse problem. Let us 
illustrate this property in a schematic example. 

We start from the zero potential with a transformation which introduces poles 
at the origin, v = diag(0, 0) — > v = diag(l, 1). In each channel we realize the usual 
(i.e. 1-channel) SUSY transformation with transformation functions sinh(Kir) and 
sinh(«; 2 r). This leads to the uncoupled superpotential 

Wd(r) = diag [m coth(fcir), k 2 coth(K 2 r)] (74) 
and the potential (see (11)) 

Vd(r) — 2diag [«; 2 cosech 2 (Kir), K 2 cosech 2 (K 2 r)] , (75) 

with the Jost solution 

, ,_ ikr( k + iK\ coth(Kir) \f k + im \ . . 

/dlK ' rj_e ^ k + iK 2 COth(K 2 r) J\ k + iK2 J [ ' 

and the Jost matrix 

F d (k) = (w doo - ikl 2 y x , Wdoo = lim w d (r) = diag (k 1 , k 2 ) ■ (77) 

r—>oo 

As coupling transformation we choose transformation function (15) with matrices C 
and D given by (37). The explicit expression for u c (r) coincides with (39) where 

= K + Kl coM«ir) e _ = . + K2 coth( K2 r) ^ r 

K + Kl K + K 2 

The parameter x from (37) should be chosen in order to avoid any singularity in the 
transformed potential. As can easily be seen from the analysis of det u c , it is sufficient 
to choose x large enough. The asymptotic behaviour of the superpotential is given by 
(40) or (41). 

The Jost matrix F c (k) may be found from (30). Its explicit expression is rather 
involved and we omit it. More important is its determinant (32), the expression of 
which is extremely simple, 

P -4- K 2 

det FJk) = ±- . (79) 
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Figure 2. Exactly solvable s — s potential V c = V c 
with parameters ki = 1.5, K2 = 1, q = 0.4, k — 
4.14286, x = 15. 



From here we find the location of the bound state at kb — in and the virtual state at 

The chain of two SUSY transformations with parameters Kip = and K, q, x 
described above leads to the mixing angle (70). The corresponding potential (q = 0.4) 
is shown in figure 2. The factorization constant k is fixed from (72). As a result, the 
mixing angle takes the form 

t 2 (k) 2fc 2 /«iK2 tana 

an e w _ Kfnl sec 2 a + k 2 (nl + k\) ' ^ ' 

Parameters k± and K2 are related with 1-channel transformations and allow us to fit the 
scattering lengths. The mixing angle t(k) depends on parameter a, which allows one 
to fit its experimental behaviour at low energies. The mixing angle at large energies 
tends to a constant value, tan 2e(fc — + 00) = —2k\K2 tan a/{n\ + «§) which can also be 
fitted using corresponding experimental data (if available) . Figure 3 shows the phase 
shifts and mixing angle for two coupled s — s potentials. 

The phase shifts of the diagonal potential Vd are shown as dotted lines in figure 
3(a). The phase shifts of the transformed potential V c are shown as dashed (q = 0.4) 
and solid (q = 1.2) lines respectively. One can see that the slopes of these curves 
coincide at the origin. The mixing angles of the transformed potential are plotted in 
figure 3(b). 

According to (79) this potential has one bound state at the factorization energy 
Eg = —k 2 . Note that the normalization constant of the bound state wave function is 
determined by parameter q as follows from (33). 

6. 3. Coupled s — p partial waves 

In this section we consider the simplest s — p coupled potential. This potential is 
characterized by I = v = diag(l, 0). The coupling transformation acts as follows: 

I = diag(Zi, l 2 ) — > I = diag(^ 2 , h) , v = diag^, v%) — > v = diag(z/i— 1, v 2 -\) .(81) 
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Figure 3. The scattering matrix for the coupled s — s potential, (a) The 
eigenphases are shown as dotted lines for Vd and as dashed and solid lines 
for V c - The parameters are: solid lines - Ki = 1.5, K2 = 1, q = 0.4, 
K — 4.14286, x — 15, dashed lines - Ki = 1.5, K2 = 1, q = 1.2, k = 13.6667, 
x = 15. (b) Mixing angle e. 



20 



V 



r-22 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 4. Exactly solvable s — p potential V" c = 
14 — 1(1 + l)r -2 with parameters Ko = 1-5, ki = 1.75, 
k = 3.53, q = 1, x = 1. 



Therefore the initial diagonal potential should have / = diag(0, 1) 
properties are satisfied for the initial potential of the form 



V d {r) = diag 



In VKlsinh (nor) , sinh (Kir)], 2r 

ar z 



(2,1). These 



(82) 



where the potential in the s-channel is obtained from the zero potential after two 
consecutive SUSY transformations with k and k,\ as factorization constants. The 
potential in the p-channel is just the centrifugal term. The Jost solution in the s- 
channel is expressed in terms of the Wronskian of factorization solutions sinh(/Cjr), 
.7=0,1, 



fd(k,r) = diag 



W [sinh (Kor) , sinh (kit) 



(k + iKo)(k + in\)W [sinh (n^r) , sinh (/tir - )] ' 



(i + kr)e 
kr 



ikr 



(83) 
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The uncoupled Jost matrix reads 

F d (fc)=diag - [{k + iK )(k + i^)}' 1 ,ik^ . (84) 

The next step is to apply the coupling transformation with the transformation 
function (15) where the Jost solution is replaced by (83). An example of potential 
curves is shown in figure 4. The Jost matrix (84) is transformed according to (30) and 
the scattering matrix is given by (46). 

Since the transformed eigenphase shifts coincide with the initial phase shifts, 
we may fit the phase shifts for the s and p waves separately before the coupling 
transformation. In our example the phase shifts read 

k k 

S c -i{k) = 0, (5 c . 2 (fc) = 7r — arctan arctan — . (85) 

Parameters Ko and k\ allow one to fit the s-channcl phase shifts. Parameter k may be 
used to fit the slope of the mixing angle (54) at zero energy. If necessary, one may use 
arbitrary chains of 1-channel transformations to get the best fit of the phase shifts. 

6.4- Coupled s — d partial waves 

The simplest s — d coupled potential is characterized by I = D = diag(2, 0). Therefore 
the initial diagonal potential should have I = diag(0, 2) and v = (3, 1). Moreover, as 
we established in section 5, 5d ; 2(0) — <5d ; i(0) = (n + l/2)ir which for n = leads to the 
following initial phase shifts dd-2(0) = k/2 and 8d-i{0) = 0. 
We start with the initial s-wave potential 

cosh 2 (K r) 

having a zero energy virtual state [3] which follows from its Jost function 

F (k) = -^—. (87) 
k + ikq 

Note that this potential and the solutions of the corresponding Schrodinger equation 
may be obtained by a SUSY transformation. This is a regular potential. To be able 
to apply the coupling transformation, we increase its singularity at the origin using 
three SUSY transformations with the transformation functions 

w(/tj, r) = Ki smh(Kir) + k cosher) tanh(/-t r) , 2 = 1,2,3. (88) 

The potential and the Jost function in the s-channel after these transformations read 

-2k 2 

Vd;i(r) = ° - -2(lnW[u(K 1 ,r),u(K 2 ,r),u(K 3 ,r)])" , (89) 

cosh (Kor) 

(fc) = (k + iKo)(k + i Kl )(k + iK 2 )(k + iK 3 ) ' (90) 

The potential in the d-channel 

6 6(3 + 6n 4 x + 64r 2 + 4^r 3 + njr 4 ) 
Vd ^ r) =^- 21UV (K4 ' r) = rH3 + 3 K4 r + K lr^ ' (91) 



Mr) = . , l2 ° (86) 



v(K 4 ,r)=c- K * r l + 3K 4 r+- ^ , (92) 
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Figure 5. Exactly solvable s — d potential V c = V c — 1(1 + l)r~ 2 with 
parameters no = 1, Ki — 1.5, K2 = 1.75, K3 = 2, K4 = 3 q = 1, x = 15, 
K — 5.53. 



5 10 




5 fc 10 



Figure 6. The scattering matrix for the coupled s — d potential. The phase 
shifts 5 c; i(fe) and 5 c -2(k) are plotted by solid lines. The mixing angle e(k) 
is plotted by the dashed line. The corresponding parameters are Ko = 1, 
Ki = 1.5, K2 = 1.75, «3 = 2, /t4 = 3 q = 1, a; = 15, ft = 5.53. 



is obtained from the centrifugal term 6/r 2 by the SUSY transformation with v(r, K4) 
as the transformation function, which decreases the singularity of the potential at the 
origin. 

The Jost solution in the s-channel is expressed in terms of the Wronskian of 
factorization solutions 



f d Jk r) = W M K i. r )> M ( K 2, r )' u (^3, ?•), /o(fc, r)} ^ 

= e , fcr fc + motanh( Ko r)_ 
fc + i«o 

The Jost solution in the <i-channel is 

.W^[w(K4,r),ft, 2 (fcr)] . . ifer / 3 \ 
/d;2(fc, r) = z— r-i — — ; ^ , h 2 (kr) = c lKr 1 + 3tn 4 r - . (95) 
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The uncoupled Jost matrix reads 



F d (k) = diag 




(96) 



(k + ino)(k + iK\){k + in 2 ){k + in 3 ) ' k 2 



which produces the eigenphase shifts 

dd:i(k) = — - } arctan — , 

/ * ' K. 



Sd-2(k) = arctan — . 

K 4 



(97) 



Next we apply the coupling transformation with the transformation function (15) 
where the Jost solution fd — diag(/d ; i, fd;2) is combined from (93) and (95). An 
example of potential curves thus obtained is shown in figure 5. The Jost matrix (96) 
is transformed according to (30) and the scattering matrix is given by (46). 

The corresponding phase shifts and mixing angle are plotted in figure 6. The 
mixing parameter e(k) is determined by (55) which, in the current case, reduces to 



We were not able to find simple expressions for the eigenphase shifts in this case. One 
can see that the mixing parameter satisfies the effective range expansion (56) (see (98) 
and figure 6). Unfortunately, this is not the case for the phase shifts (see figure 6). 

7. Conclusion 

In the present paper, a careful analysis of SUSY transformations between uncoupled 
potentials with equal thresholds but arbitrary partial waves in N channels and coupled 
ones is given. In particular, we formulated conditions imposed on the transformation 
function to get a nontrivially coupled scattering matrix. A family of iso-phase 
potentials generated by a coupling SUSY transformation has been obtained. This 
family is parameterized by a real symmetric M x M, < M < N, non-singular matrix 
Y . The analysis of the zeros of the Jost-matrix determinant for these potentials 
showed that the SUSY transformation creates a new M fold degenerate bound state 
energy Eb = —k 2 and an N — M fold degenerate virtual state energy E v = —k 2 . 

In the most important practical case, the two-channel case, we analyzed the 
behaviour of the supcrpotcntial and potential at large distances in details. We have 
found an unusual effect, i.e. a modification of the long-range behaviour of the potential 
under a coupling SUSY transformation which consists in the exchange of attribution 
of the partial waves between the channels. Moreover, we analyzed in details the 
transformed scattering matrix, phase shifts and mixing angle. The analysis of the 
phase shifts and mixing angle demonstrated how scattering properties change after a 
SUSY transformation. 

As an illustration of our approach, several simple examples have been presented. 
First, to emphasize the difference between couplings in the potential, Jost and 
scattering matrices, we presented examples of a trivially coupled scattering matrix 
corresponding to non trivially coupled potential and Jost matrices. These examples 
answer the general questions raised in the introduction: situations may exist where 
a non trivially coupled potential leads to a trivially coupled S'-matrix, with either a 
trivially or non trivially coupled Jost matrix. Thus, the requirement that an S'-matrix 
is non trivially coupled is more restrictive than the similar requirement for a potential 




(98) 



Coupling between scattering channels with SUSY transformations 



19 



matrix. Moreover, even the requirement that the Jost matrix is non trivially coupled 
is less restrictive than the corresponding requirement for the S'-matrix. 

Afterwards, a non trivial coupling was introduced in the s — s, s — p and s — d 
channels. In both s — s and s — p examples, we have shown how to fit the low-energy 
behaviour of the phase shifts and mixing angle using parameters of the transformation. 
We note that the s — p example contains all essential ingredients of a convenient 
inversion scheme since the s — p coupling transformation conserves the initial phase 
shifts, which allows separating the inversion of phase shifts from that of coupling 
parameters. In the s — d case, to satisfy the effective range expansion for the mixing 
parameter, we used an initial potential with a zero energy virtual state. Nevertheless, 
the obtained phase shifts of the coupled s — d potential do not satisfy the correct 
effective range expansion. Moreover, the presence of the zero energy virtual state 
strongly restricts possible applications of our method to the inversion in this case. 
We believe that a chain of two coupling transformations may allow avoiding this 
inconvenience. 
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